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Abstract 

We calculate the cross section for the 77 — > J/tpJ/ip process. Two mechanisms are consid- 
ered: box (two- loop) diagrams of the order of 0(a 2 m a 2 ) and two-gluon exchange of the order of 
0(a^ m af). The first mechanism is calculated in the heavy-quark non-relativistic approximation 
while the second case we also include the effects of quantum motion of quarks in the bound state. 
The box contribution dominates at energies close to the threshold (W < 15 GeV) while the two- 
gluon mechanism takes over at W > 15 GeV. Including the bound-state wave function effects for 
the two-gluon exchange mechanism gives a cross section 0.1 - 0.4 pb, substantially smaller than 
that in the non-relativistic limit (0.4 - 1.6 pb). We also find a strong infrared sensitivity which 
manifests itself in a rather strong dependence on the mass for the t-channel gluons. The elementary 
cross section is then used in the Equivalent Photon Approximation (EPA) in the impact param- 
eter space to calculate the cross section for 208 Pb + 208 Pb — > 208 Pb + J/ipJ/ijj + 208 Pb reaction. 
Distributions in rapidity of the J/i/jJ/i/j pair and invariant mass of the pair are shown. 

PACS numbers: 12.38.Bx, 13.85.Ni, 14.40.Pq 
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I. INTRODUCTION 



The two-photon collisions is a natural place to study vector-meson pair production. Pro- 
duction of light vector mesons is rather of nonperturbative nature and therefore subjected 
to more phenomenological studies. Heavy vector meson production (J/ip or T) is particu- 
larly interesting as here the pQCD degrees of freedom can be applied and hopefully reliable 
predictions can be obtained. It was realized quite early that in the high-energy limit the two- 
gluon exchange is the dominant reaction mechanism of two heavy vector meson production. 
The relevant amplitude was calculated first in the heavy-quark non-relativistic approxima- 
tion in [l|, and first estimates of the cross section were presented there. The impact factors 
of [l| were used next to estimate BFKL effects e.g. in (3, 0|. Some other estimates based 
on parametrizing dipole-dipole interaction were presented in j^]. Perhaps surprisingly, there 
is actually a substantial spread in the predictions shown in those works, they do not seem 
consistent one with each other and differ by almost two orders of magnitude. For a number 
of reasons, it seems that these previous evaluations are not very realistic in the intermediate 
energy range relevant to studies of the 77 sub-process in nucleus-nucleus collisions. 

Firstly, the process considered in the mentioned papers is not the lowest order of strong 
coupling constant. Formally lower-order processes of the box type were studied in jBf. 
In these calculation heavy-quark approximation was used in order to calculate the two- 
loop amplitudes. At high energies, these mechanisms correspond to the quark- ant iquark 
exchanges in the crossed channels and will die out with increasing 77 energies. The numerical 
predictions of 0] were however shown only at high (linear collider) energies (a/s > 400 GeV) 
where this mechanism is definitely not the dominant one. 

Secondly, previous calculations of the two-gluon exchange mechanism were restricted to 
the extreme non-relativistic limit, and neglect the motion of heavy quarks in the bound 
state. Here, we shall calculate the two-gluon exchange contribution using for the first time 
for this process relativistic non-forward impact factors. Furthermore, by introducing a gluon 
mass, we investigate the infrared sensitivity of the amplitude, which turns out to be rather 
strong. 

Therefore in the present work we wish to discuss both mechanisms simultaneously and 
discuss and identify the region of their dominance. 

Finally we shall present realistic predictions for heavy ion collisions where the fluxes of 
quasi-real photons are very large. Here one may be able to study such processes in the near 
future. 



II. FORMALISM 



In the present approach we include processes shown in Figs, [fl [21 El Now we shall discuss 
each of the mechanisms one by one. 



A. Box diagrams 

At the leading order, (D(a 2 a 2 ), the sub-process 7 + 7 — > J/ij) + J/ij) is represented by 
the 20 "box" diagrams shown in Fig{T]The calculations are straightforward and follow the 
standard QCD rules and non-relativistic bound state formalism described in detail in Refs. 
0. As usual, the production amplitudes contain spin and color projection operators that 



2 




FIG. 1: The "box" diagrams for 77 —> J/vpJ/ip included in the present paper in the heavy-quark 
approximation . 
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FIG. 2: Graphical representation of the two-gluon exchange mechanism for 77 — > J/ipJ/i/j. 

guarantee the proper quantum numbers of the final state mesons. The formation of the 
cc bound states is determined by the J/ip wave function at the origin of coordinate space, 
|^(0)| 2 - its value is known from the leptonic decay width: (^(O)] 2 = |7?.(0)| 2 /(47r) = 0.08 
GeV 3 [2J. Only the color-singlet channels are taken into consideration in the present study. 

Our calculation is identical to that of Ref . [5] , with the only exception that we consider 
individual helicity amplitudes rather than the entire spin- averaged matrix element. The 
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FIG. 3: Graphical representation of the meson exchange diagrams for 77 — > J/i/jJ/i/j. 

polarization vectors of the initial photons and outgoing mesons are presented in our approach 
as explicit 4- vectors. For a particle moving along the z-axis, = (0, 0, \p\, E), the helicity 
eigenstate vectors have the form 

e"(±l) = (±1, i, 0, 0)/V2, e"(0) = (0, 0, E, \p\)/m. 

The evaluation of the Feynman diagrams has been done using the algebraic manipulation 
system FORM We have checked that after performing numerical summation over all 
P_ossible polarization states we arrive at the same result as the one given by Eq.(4) in Ref. 

In practical calculation we take leading-order formula for a s which is evaluated at fi 2 . = 
Am 2 r . 



B. Two-gluon exchange 

At sufficiently high energies, the box diagram contributions, which contain the fermion- 
antifermion exchange in the crossed channels die out, and the cross section will be dominated 
by diffractive mechanisms. The Feynman diagrams for the diffractive 77 — > VV amplitude 
is depicted in FigfSJ Although it is formally of higher order in as than the box mechanism, 
the crossed channel gluon exchanges do not die out with energy. 



1. Relativistic approach 

The altogether 16 diagrams of the type shown in FigJH can be lead to the amplitude, 
which can be cast into the impact-factor representation: 

Here q is the transverse momentum transfer, t ~ — q 2 , and \Lq is a gluon mass parameter. 
Notice, that the amplitude is finite at ji G — > 0, because the impact factors J vanish for 
k, — > ±q/2. Graphical representation of the impact factors is shown in FigJH 
The amplitude is normalized such, that differential cross section is given by 

^iSK 7 ^^^ • (2 - 2) 
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At small t, within the diffraction cone, the cross section is dominated by the s-channel 
helicity conserving amplitude. In this case, the explicit form of the impact factor is 



J{l\ -»• K; q) = $\,t V N* - 1 a s e Q V^a em J 



(L(l2k Mz,k)I(T,T), (2.3) 



z(l - z)(2tt) 3 



where N c = 3 is the number of colors, cq is the charge of the heavy quark, and ipv{ z i is the 
light-cone wave function of the vector meson. It depends on the light-cone-plus momentum 
fraction z carried by the quarks as well as on the quark transverse momentum k. The 
spinorial structure of the V — > QQ vertex has been chosen such that it describes the S-wave 
bound state of quarks [9], so that 

I(T, T) = m 2 Q $ 2 + [z 2 + (1 - zf] (fc#0 + M ™ Q 2m [fc 2 $ 2 - (2z - l) 2 (fe$i)] . (2.4) 

Here #i,$2 are shorthand notations for the momentum structures, corresponding to the 
four relevant Feynman diagrams: 

1 1 1 1 

$9 = 1 h 

* / 1 i \9 , 9 /I \9 , 9 / J i /r-i\9 , 9 



(I + K )a + m a (Z_ K )2 + m ^ (Z + q /2 )2 + m ^ (i _ q/ 2 )2 + m 2 Q 

1 (Z + ^ + m 2 , (Z-K)2 + m 2 + (i + q /2)2 + m 2 Q + (Z-q/2)^ + mJ' 
and we used 

/ 1 \ k 2 + rrin 

l = k+ (*-2)«' M2 =m^- (2 - 6) 

The electronic decay width T(V — > e + e~) of the vector meson constrains the light-cone wave 
function. Explicitly, the decay constant gy is 

/dzd 2 k 2 
«(1 - z)(27rf 3 

it is related to the decay width through 

Ana 2 „ 

r(y^e + c -) = _^«.^. (2.8) 

In the present approach we assume real photons and therefore only transverse photon and 
vector- meson polarizations are taken into account. This is sufficiently good approximation 
for heavy-ion peripheral collisions where the nucleus charge form factor selects quasi-real 
photons. 

The running scale of strong coupling constant for the evaluation of the two-gluon exchange 
cross section is taken as: 

2 r 2 jr 2 i 2 l 

\i = maxj/t , k + ■ 

(2.9) 
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For the radial wave function we use a Gaussian parametrization 

1 / k 2 4~ rn? \ 
1>v{z, k) = N ■ exp[-ap 2 ] , p 2 = - { ^ _ ^ ~ 4m 2 ,) . (2.10) 

Where the constants N, a are determined from the decay width (12.7112.81) and the normaliza- 
tion condition 

^7(^^) M2|M2 ' fe)|2=1 ' (2 - u) 

For the case of J/ip of interest here, we use mg = m c = 1.5 GeV and a = 1.37 GeV -2 . In 
this way we get a good description of the jp — > J/ipp experimental data (lo| . 



C. Non-relativistic limit 



In previous calculations in the literature the extreme non-relativistic approximation ap- 
propriate for the weakly bound state of heavy quarks has been adopted. In this case any 
quantum motion in the bound state is neglected, and one uses effectively: 

ipv(z,k) = C-5(z- ^)5 (2) (fc). (2.12) 

Notice that this means that quark and antiquark momenta in the boosted bound state are 
collinear. The constant C results as 



NJ4 V 7T w e Q a cm N c M v \l 71 1 ' 

Furthermore, for the weakly bound state My = 2rriQ. In the impact factor, we can now 
approximate I(T,T) = mq$ 2 , and the structure $2 simplifies to (hereafter NR means non- 
relativistic): 



k 2 + M 2 /4 q 2 + M 2 



* 2 NR = H ... , —\- (2-11) 

The impact factor thus becomes: 



J^ lx -> V T ;n,q) = , v ^KE± a s f . 



D. Meson exchange processes 

The J/?/> meson decays into a photon and pseudoscalar meson [7]. The corresponding 
partial decay width can be used to calculate corresponding coupling constant g 2 MJ /^ oc 

^J/ip-^-yM(0-)/Pdec- 
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FIG. 4: Graphical representation of the two-gluon impact factor for the 7 — > J/tp transition. 

The amplitude for the 77 — > J/ipJ/ip due to pseudoscalar meson (M) exchange shown 
in Figj3]can be written as: 



FjMJ/if> —F'yMJ/i) (t) 



— I 

l M 

+ g-yMj/^^e^ (Ai)e^ (A 4 )&i jai p 2A 

— 1 

FjMJ/ip (u) —F^Mj/ip (u) 

u m M 

g^Mj/^ a ^e, 2 (X 2 )e: 2 (\ 3 )k 2 , a2 p 1A . (2.16) 



Above we add t and u meson exchange amplitudes. Compared to the J/tp — > 7M decay in 
the 77 — > J/ijjJ/i/j process the exchanged mesons (M = 7r,r],r] c ) are off mass shell. This 
requires to introduce corresponding form factors (F in (12 . 1 6[) ). In the following they are 
parametrized in the exponential form as: 

F(k*) = exp (^^) , (2.17) 

where k is four-momentum of the exchanged meson. We take A = 1 GeV to estimate the 
meson exchange contributions. We have found that the cross section for the exchange of 7T°, 
r) and rj c is very small, several orders of magnitude smaller than the cross section for the box 
and two-gluon exchange mechanisms. This is very different than for the production of pairs 
of light mesons in photon-photon scattering where the exchange of t or u channel mesons 
lead to relatively big cross section of the order of 1-10 nb. 

The contributions of 7r and rj meson exchange are small because corresponding coupling 
constants J/ip- meson-photon are small compared to the case of p or u decays as can be 
deduced from the partial decay widths. The contribution of rj c exchange is suppressed by 
large rj c mass in the propagator. 



E. Equivalent Photon Approximation for AA — > AAJ/ipJ/tfj 
The total cross section can be calculated by the convolution: 
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a (AA -> Jf^Jf^AA- s AA ) = [ d Xl dx 2 a (77 #M ^2^) ^77 ^ij^,b) 

7 dx\dx 2 

The effective Weizsacker- Williams photon fluxes can be calculated from: 
rfn 77 (x 1; x 2 , b)= f d 2 b x d 2 b 2 5 (2 \b -bt + b 2 )Sl bB {b) -^\E(x u b Y )\ 2 -^l\E(x 2 , b 2 )\ 2 

J ' IT Xi 7T X 2 

(2.19) 

The presence of the absorption factor S 2 hs (b) assures that we consider only peripheral 
collisions, when the nuclei do not undergo nuclear breakup. Following 11] this can be taken 
into account as: 

5 a 2 bs (b) =9(b- 2R A ) = 6 (lb, -b 2 \- 2R A ) . (2.20) 

Thus in the present case, we concentrate on processes with final nuclei in the ground state. 
The electric fields can be expressed through the charge form factor of the nucleus: 

E (X, b) = ZV^Z J -0-f- lQb g2 + l2 M 2/ em {<? + ^ M\) . (2.21) 

Different forms of charge form factors are used in the literature. Here we use charge 
form factor calculated as a Fourier transform of the charge distribution. It was shown in 
Ref. 12, 13 to be crucial to take the realistic form factors. 



III. RESULTS 

A. 77 — > J/tpJ/ip reaction 

Let us start from presenting our relativistic impact factors. In Figj5]we show the forward 
(q 2 = 0) impact factor as a function of gluon transverse momentum squared. 

One can observe a quick rise of the impact factor till gluon momentum transfer k 2 w 
3 GeV 2 . The impact factor obtained in the present relativistic approach (solid line) is 
much smaller than that obtained in the non-relativistic approach (dashed line) used in the 
literature. This will lead to corresponding differences in the amplitude of the 77 — > J/i/jJ/i/j 
process and related cross section. 

Let us illustrate the situation for the amplitude. In Figj6] we show rather integrands: 
dA/dn (left panel) and dA/dtfi (right panel) instead of the whole amplitude. Here is 
azimuthal angle between two-dimensional vectors k and q . This figure demonstrates also 
that it is relatively easy to get numerical convergence of the integration in calculating the 
amplitude of the process, the decrease in k is rather fast, the oscillation in is relatively 
shallow. 

In FigJ7]we show the dependence of the cross section on two- momentum transfer squared 
at energy y/s ^> 2mj/^. These two-gluon exchange distributions scale at larger energies. 
At smaller energies kinematical limitations due to energy-momentum conservation as well 
as interference of t and u diagrams due to identity of particles (two identical J ftp's) have 
to be included which destroys the scaling behavior. Fortunately both effects happen when 
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FIG. 5: Forward impact factor for 7 — )■ J/ip as a function of gluon momentum transfer squared. 
The solid line is our new result within the relativistic approach. Here ^j/^(K 2 )S\ 1 \ a = Ji^lXx — > 
V\ 3 ) K, q = 0) (see Eq. (|2.1|) ). For comparison we also show the impact factor in the non-relativistic 
approach (dashed line). 




k (GeV) cp (deg) 



FIG. 6: Integrands of the amplitude as a function of k for = (left panel) and as a function of 
(f) for q% = 5 GeV (right panel). The solid line is our new result within the relativistic approach. 
Form comparison we show results for the non-relativistic approach (dashed line). 

W < 10 GeV, i.e. when the contribution of the two-gluon exchange is much smaller than 
that for the box mechanism as will be shown below. The power-like distribution shown in 
the figure strongly deviates from exponential dependence which was assumed for simplicity 
in previous calculations (see e.g. (1,0]). 

A brief discussion of the infrared sensitivity of the cross section is in order. Certainly, the 
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cross section is infrared safe - the impact factors vanish, when the transverse momenta of t- 
channel gluons go to zero. However this does not mean that there is no numerical sensitivity 
to the domain of small gluon transverse momenta, where the use of perturbation theory 
is not warranted. To investigate this sensitivity, we introduced a mass for the t-channel 
gluons. Notice that this mass does not affect the gauge invariance driven cancellation of 
impact factors which guarantee finiteness of the cross section. 
Numerically, the infrared sensitivity is quite strong: 

ff(77 ^ j / ^j M/iG = 0, 7G eV) 

<7( T 7^ JM/M/^ = 0GeV) ~ ' ' 1 ' ) 

Here, the value \i& = 0.7 GeV corresponds to a vacuum correlation length of gluons of 
R c ~ 0.29 fm. It is motivated by a number of phenomenological approaches to high-energy 



scattering |14|, |15j|. An analysis of gluon field-strength correlators on the lattice |16| also 
suggests a correlation length in this ballpark. Our observation of a surprisingly strong 
infrared sensitivity of the 77 — > J/ipJ/ip cross section agrees with Ref. 17J, which uses the 
non-relativistic limit of the amplitude. 




FIG. 7: Dependence of the 77 — > J/ipJ/ifj cross section on momentum transfer squared. The solid 
line represents relativistic approach and the dashed line the non-relativistic approach. 

The total cross section for the 77 — > J / ipJ / ip reaction is shown in FigJHJ We show separate 
contributions of the box diagrams (dashed line) as well as of the two-gluon exchange (dotted 
line). The second contribution scales at larger energy. Any higher-order interaction (t- 
channel ladder exchange) would increase the Born cross section. Unfortunately the existing 
calculation in the literature are limited to leading-order BFKL type of calculation 0, y| 
which is known at present to strongly overestimate the cross section for other reactions. We 
leave a realistic evaluation of the increase of the cross section for future studies. 

The box and two-gluon exchange cross sections have quite different energy dependence. 
As a consequence they are of similar size only in a very limited range of energy. Therefore 
in practice the interference effect can be neglected. 

The suppression of the relativistic amplitude in comparison to its non-relativistic limit 
can be traced back to the suppression from the off-shell quark propagators in Eq. fl2.5p . 
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FIG. 8: Energy dependence of the 77 —> J /if) J /if) cross section. The solid line is for the the 
box contribution (dashed line) while the dashed and dash-dotted lines for two-gluon exchange 
contribution in non-relativistic and relativistic approaches, respectively. 



B. AA -> AJ/ipJ/ipA reaction 

Now we pass to the calculation of the cross section for nuclear collisions. We follow 
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cc 



18 



Or 7T7T 



19 . In those 



similar analyses for exclusive production of p°p° [l2j, 
approaches a calculation of the distribution in rapidity of the pair or invariant mass of the 
pair is particularly easy. In Figj9] we show distribution in rapidity of the pair for Wnn 
= 2.76 TeV. Both contributions (box and two-gluon exchange) are shown separately. The 
phase space integrated nuclear cross section is a = 0.1 ph. which is of similar size as the 
cross section for production of DD meson pairs [20J]. The contribution of the two-gluon 
exchange is an order of magnitude smaller than that for the box mechanisms. No previous 
evaluation of the nuclear cross section included the dominant box contribution. In FigJTUIwe 
show distribution in J ftp J /if) invariant mass. This distribution drops quickly as a function 
of the two-meson invariant mass, which reflects the fall-off of the effective photon-photon 
luminosity at large 77-energies. The two-gluon exchange contribution shows up only at 
Mj/ipj/f > 15 GeV. The situation can improve slightly if the rise of the 77 cross section due 
to BFKL-type effects were included. We do not expect, however, that a realistic calculation 
including a BFKL-type ladder interaction between two cc pairs would change our predictions 
significantly for invariant J /if) J /if) masses smaller than say 50 GeV. 

Recently we have calculated the cross section for exclusive production of single J /if) 
meson [2l[. In Fig{TTJwe compare the cross section for the pair production (photon-photon 
fusion) with that for the production of a single J /if) via the "photon-pomeron" mechanism. 
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FIG. 9: Rapidity distribution of the J/ipJ/if> pairs in 208 Pb + 208 Pb -^ 20S Pb + J /if) J /if + 208 Pb 
process for Wjvat = 2.76 TeV. The dashed line is for the box contribution while the dotted line 
corresponds to the two-gluon exchange. 




FIG. 10: Invariant mass J/ifJ/if pair distribution in 208 Pb + 208 Pb -> 208 Pb + J/ipJ/ip + 208 Pb 
process for Wnn = 2.76 TeV. The meaning of the curves is the same as in the previous figure. 

The latter cross section is more than five orders of magnitude larger. The single J /if} 



production was measured e.g. at RHIC [22]. Similar analysis is being performed by the 
ALICE collaboration at the LHC [23]. A present statistics is a few hundreds of J /if) [24] . 
This means that with present statistics the number of measured pairs would be less than 
one. Clearly a better statistics is necessary to see pairs of J/ip. A better situation could 



12 



be for the pp — > ppJ/ipJ/ip reaction where the effective 77 luminosity at large J/ipJ/ip 
invariant masses is much larger. 

If the nuclear cross section for 208 p5-|_ 208 p5 __j.208 p^j j^pj j^Y^Pb was much larger than 
predicted in this paper it could mean that a double scattering photon-pomeron mechanism 
discussed in Ref.[2l|] plays important role. A measurement of pair production at ALICE 
can therefore provide a quite new information. The exclusive nuclear double-scattering 
production of mesons was not discussed so far in the literature. 
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FIG. 11: Energy dependence of the cross section for the 208 P6+ 208 Pb -> 208 Pb+ J/tpJ/ip + 208 Pb 
process versus that for the 208 Pb + 208 Pb — > 208 Pb + <//?/> + 208 Pb process. The latter cross section 
is taken from [211 ] . 



IV. CONCLUSIONS 

In the present paper we have calculated the total cross section and angular distribu- 
tions for the 77 — > J /ip J ftp process. For the first time we have included both formally 
two-loop contribution, called here "box contribution", included so far in only one previous 
study, as well as two-gluon exchange contribution, formally higher order, discussed already 
in the literature. We have clarified a disagreement of results of different calculations in 
the literature. In our approach the two-gluon exchange contribution has been calculated 
in a relativistic approach with the wave function adjusted to reproduce the HERA data 
for 7£> — > J /ipp reaction. For the first time we have calculated relativistic off- forward two- 
gluon exchange amplitude. We have quantified a difference between the relativistic and 
non-relativistic approach used so far in the literature. The relativistic approach with "real- 
istic" quark- ant iquark J/i/j wave function gives smaller cross sections than those obtained in 
non-relativistic approximation. In addition, we have also estimated a contribution of pseu- 
doscalar (tc°, rj, r] c ) meson exchange mechanism. The latter ones turned out to be practically 
negligible. 
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The elementary cross section has been used to calculate cross section for ultra-peripheral 
ultra-relativistic heavy ion collisions (AA — » AJ/if)J/ipA) at the LHC In this context we 
have used equivalent photon approximation in the impact parameter representation and 
realistic charge form factor of the 208 Pb nucleus. The nuclear cross section has been estimated 
and distribution in the J /if) J /if) rapidity and invariant mass have been presented. The cross 
section for double J /if) production was compared with that for single J /if} production from 
(2l| . At the LHC energies the double J /if? production cross section is smaller by more than 
five orders of magnitude compared to that for single J /if? production. A measurement of 
the J /if) J /if> pair production may be then rather difficult as one has to take into account 
in addition a product of respective branching fractions for decays into lepton-antilepton 
pairs. On the other hand much bigger cross section could be a signal of double scattering 
"sequential" production of two single J/ip's. 

In the present analysis we have concentrated on the photon-photon mechanism of the 
J /if) J /if> pair production. Evidently similar mechanisms exist for gluon-gluon (sub)processes, 
leading to the inclusive production of J /if) J /if) pairs in proton-proton collisions. Such a 
production has been observed recently by the LHCb collaboration 25[. Only the box-type 



diagrams were included in corresponding calculations [26l . |27| . The two-gluon exchange 



mechanism considered here may in principle compete with the double-parton production 



mechanism of the J/if?J/if> pairs discussed recently in [26l . |27 . 

After we have performed the calculation presented here we have found a new preprint [28| 
where the production of two vector mesons is discussed. Some additional photon-exchange 
processes are considered. In the 77 channel, the authors do not discuss the dominant box- 
mechanism. Also, nuclear charge form factors which are important at large 77 energies 



12l . Il3| are neglected. 
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